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Abstract. 

Over the last years the exciting developments in the field of ultracold atoms confined 
in optical lattices have led to numerous theoretical proposals devoted to the quantum 
simulation of problems e.g. known from condensed matter physics. Many of those 
ideas demand for experimental environments with non-cubic lattice geometries. In this 
paper we report on the implementation of a versatile three-beam lattice allowing for 
the generation of triangular as well as hexagonal optical lattices. As an important step 
the superfluid-Mott insulator (SF-MI) quantum phase transition has been observed 
and investigated in detail in this lattice geometry for the first time. In addition to 
this we study the physics of spinor Bose-Einstein condensates (BEC) in the presence 
of the triangular optical lattice potential, especially spin changing dynamics across 
the SF-MI transition. Our results suggest that below the SF-MI phase transition, 
a well-established mean-field model describes the observed data when renormalizing 
the spin-dependent interaction. Interestingly this opens new perspectives for a lattice 
driven tuning of a spin dynamics resonance occurring through the interplay of quadratic 
Zeeman effect and spin-dependent interaction. We finally discuss further lattice 
configurations which can be realized with our setup. 
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1. Introduction 

The physics of quantum degenerate gases in optical lattices has rapidly grown to one 
of the most accounted fields of atomic physics over the last couple of years [I]. This 
new experimental environment opens the way to study many-body physics from so 
called mean-field dominated physics to strongly correlated regimes and allows for a 
vast amount of interesting experiments ranging from Josephson junctions [2J, quantum 
registers [3] and ultracold chemistry [I] to the implementation of solid state models 
[5] and experiments concerning spin oscillations in deep optical lattices [6] to name 
only a few examples. However, up to date nearly all experiments with ultracold atoms 
have been performed in simple cubic lattices owing to the ease of their experimental 
implementation as compared to other lattice geometries [7]. It is only recently that 
the experimental and theoretical interest in more complex lattice configurations such as 
arrays of fully controllable double wells [3 [9] or artificially generated gauge potentials 
[TOl [TTj appears. In particular it has been proposed p21 031 EH EEJ HE] that cold atoms 
loaded into triangular or hexagonal optical lattices exhibit very rich new phases showing 
e.g. supersolidity [HI [12], quantum stripe ordered states [15], exotic superconducting 
states [16] or graphene-like physics [UJ, [18]. In solid state physics frustrated quantum 
magnets on triangular or trimerized lattices attract the highest amount of interest due to 
their intrinsic link to high temperature superconductivity (see e.g. [19] and references 
therein). Recently it was proposed how frustrated quantum magnets can be realized 
within a triangular optical lattice |20j. 

In the context of magnetic interaction in ultracold quantum gases, spin mixtures are 
currently receiving rapidly growing experimental and theoretical interest. They combine 
the unprecedented control achieved in single component Bose-Einstein condensates 
(BEC) with intrinsic degrees of freedom. The spin-dependent coupling connects 
quantum gas physics to fundamental magnetic phenomena. Earlier work analyzed the 
ground states and dynamical evolution of 23 Na [2D [22] and 8?Rb [23 123 123 123 [23 123129] 
spinor condensates in various ground state hyperfine manifolds. Fascinating phenomena 
like spontaneous symmetry breaking [30], dynamical instabilities leading to pattern 
formation [31j [32] as well as dipolar effects [33J have been observed in harmonically 
trapped spinor condensates. However experiments dedicated to the investigation of 
multicomponent systems in optical lattices have been restricted to isolated 2- and 3- 
body spin changing oscillations deep in the Mott-insulating regime which are inherently 
different from the mean-field driven dynamical evolution in harmonically trapped BEC 
P, [34] . The physics of spinor condensates in the presence of periodic potentials but still 
in the superfluid regime is thus far widely unexplored. In this manuscript we report on 
a new experimental setup realizing a 87 Rb BEC in a triangular optical lattice [?J. 
We focus on a concise understanding of the Mott-insulator (MI) transition in the 
triangular lattice which is a necessary precondition for further experiments on quantum 
magnetism in weakly as well as in strongly correlated regimes. In particular we work 
in a 2D system, where the higher dimensional degrees of freedom are frozen out by a 



Ultracold quantum gases in triangular optical lattices 



3 



sufficiently tight transverse confinement mediated by a ID lattice. 

Furthermore we investigate spin dynamics of 87 Rb F — 1 spinor BEC in this new lattice 

geometry. 

2. The three-beam lattice 

We start with a description of the experimental realization of the triangular optical 
lattice. The periodic potential employed in our experiment is created by three laser 
beams which intersect in the x — y plane mutually enclosing angles of 120° with lattice 
vectors ki = k (1, 0), k 2 = k (-• v/3/2, -1/2) and k 3 = k (y/3/2, -1/2) (Figfljt). k is the 
wave vector of the laser employed for the lattice beams. Adding up the field strengths 
of the individual beams coherently 
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E 2D (r, t) = Y^ E oi e< cos(k r r -ut + fa) (1) 

i=l 

and averaging over time results in a potential that can be written as 

/3 1 \ 
V(r) = V I - + - (cos (bir) + cos (b 2 r + 12 ) + cos ((b x - b 2 )r + 23 ))J (2) 

for the case that all the same i.e. the polarization is perpendicular to the plane 

spanned by the lattice beams. For simplicity the relative phases faj will be set to zero 
in the following without loss of generality. Vo = 4Zo corresponds to the depth of a ID 
lattice created by two counterpropagating beams, where Iq is the light shift produced 
by one of the lattice beams. The corresponding reciprocal lattice vectors bj can simply 
be derived through the relation bj = e^k^kj — kfc), while the primitive direct lattice 
vectors follow according to a^ • bj = 2ii5ij. The resulting potential is plotted in FigHK 
and obeys a triangular symmetry. It is important to note, that in this configuration 
the polarization is purely linear and the potential therefore almost m^-independent for 
large laser detuning. A variation of the phases 0« of the laser beams only affects the 
lattice through a global shift in position. This very fact can be readily used to move, 
accelerate or even rotate the lattice [20] . In order to directly illustrate the effect of the 
triangular lattice structure at this point, Fig. [2^ shows absorption images which have 
been obtained after switching off all optical potentials abruptly which projects the quasi 
momentum distribution in the lattice on free momentum states. Neglecting interaction 
effects the density distribution after a sufficiently long time-of-flight (TOF) therefore 
represents the quasi momentum distribution in the lattice [35J which directly reveals 
the reciprocal lattice of the triangular lattice is shown in Figj2^,. Opposed to that Fig. 
Eb reveals the shape of the first Brillouin zone which can be obtained by first populating 
the whole first band by carefully exciting the atoms followed by an adiabatic ramp-down 
of the lattice potential thereby mapping the quasi momentum on free momentum states. 
By rotating the polarization of the three lattice beams in the x — y plane as depicted in 
Fig [lb, another lattice geometry can be realized. The potential minima are now ordered 
on a hexagonal lattice as shown in Fig [lb. Moreover the polarization at the potential 
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Figure 1: Three beams intersecting in a plane and enclosing angles of 120° pairwise 
create a triangular optical lattice if their polarization is perpendicular to the plane 

(a) . When the polarization is rotated into the plane spanned by the lattice beams 

(b) a hexagonal lattice with alternating circular polarization at the potential minima 
is obtained. A perpendicular ID lattice provides an overall three-dimensional periodic 
confinement. 

minima is perfectly circular and alternating in helicity for nearest-neighboring sites. For 
alkali atoms, such as 87 Rb the dipole force can be rewritten in a very instructive form 
if the laser detuning A is on the order of the fine-structure splitting Aps in 87 Rb 

^« = 10( 1 + ^ir)'«< < 3) 

where u is the laser frequency, nip the magnetic sub-state of the atoms and V = 0, ±1 
characterizes the laser polarization (ir, a ) respectively. This gives clear indication that 
for circular polarization the dipole potential will obey a spin dependence, which vanishes 
asymptotically for infinite detuning. Taking this into account, the hexagonal lattice will 
induce an anti-ferromagnetic ordering for atoms with nip ^ . To obtain a quantitative 
expression for the hexagonal potential it is convenient to decompose the electric field 
into its components in a circular basis e CT ±. The corresponding potential for the two 
polarizations is thus given by 
V 

V + = — (3 + 2(cos(bx ■ r — 2(f> c ) — cos(b 2 • r — C ) + cos((b! — b 2 ) ■ r — 2<p c 
8 

Vn 

V- = (3 + 2(- cos(b! ■ r - C ) + cos(b 2 • r - 2(p c ) + cos((bi - b 2 ) • r + 2 
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with the characteristic phase <p c = ir/3. 

These two expressions clearly emphasize, that the translational symmetry of the 
underlying sub-lattices is exactly the same as for the triangular lattice. The change in 
polarization only corresponds to a change of the local basis with respect to the position 
of the individual potential wells inside a primitive cell. 

By finding the maxima of the above expressions one can easily show that the primitive 
(a) (b) 
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Figure 2: (a) Absorption images of ultracold atoms released from a triangular optical 
lattice. The sample has been imaged from two different directions (b) Result of a full 2D 
band structure calculation for the triangular optical lattice for a lattice depth of 3 E-&. 
The lowest four bands are presented along lines connecting points of high symmetry as 
depicted in the inset. An experimental time-of-flight image of the Brillouin zone is also 
shown. 

cell of the hexagonal lattice contains an anti-ferromagnetic basis with one <7 + -well at 
a* = 1/3 (ax + 2a 2 ) and a cr~-well at a£ = 1/3 (2a x + a 2 ). Disregarding the nip- 
dependence as justified for \F,0), the pure intensity modulation given by the sum of the 
two terms in [5] can be written as 

^2D(r) = V (| - ~ (cos (bi • r) + cos (b 2 ■ r) + cos ((b a - b 2 ) • r))) . (6) 

For other m^-states of 87 Rb in F = 2 and a laser wavelength of 830 nm the relative 
difference in the optical dipole potential can be as large as (Ifj^^ ~ ^dip^V^iip ~ 20% 
according to Equ. [3j We emphasize that the interplay between spin-dependent contact 
interaction, frustrated inter-site dipolar interaction and a spin-dependent trapping po- 
tential where e.g. spin-dependent tunneling could be realized may give rise to a whole 
range of new physical effects and exotic ground state phases. 

3. Experimental preparation of BEC and lattice calibration 

In the following we briefly summarize the experimental sequence employed for the mea- 
surements presented here. BEC is produced by first collecting up to 10 10 87 Rb atoms 
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in a double magneto-optical trap system, which are subsequently loaded in a magnetic 
trap of hybrid D-cloverleaf type. Evaporative cooling leaves us with about 1.5 x 10 6 
ultracold atoms slightly above the critical temperature for Bose-Einstein condensation. 
Afterwards the atoms are transferred to a crossed dipole trap with typical trapping fre- 
quencies uj = 2tt x (90±3)Hz. The large detuning of the employed Nd:YAG laser with a 
wavelength of A = 1064 nm ensures a spin state independent global trapping allowing for 
the investigation of spinor physics at all lattice depths. Further evaporative cooling by 
smoothly lowering the intensity of the dipole trap laser beams produces BEC with parti- 
cle numbers ranging from 6 to 7x 10 5 without any discernible thermal fraction. Since we 
are especially interested in the physics induced by the triangular lattice we eliminate the 
influence of the third dimension by creating a stack of quasi two-dimensional disc-shaped 
condensates. This is done by ramping up a ID lattice over 150 ms to a final depth of 
Vo = 30E r , where the recoil energy is given by E r = (ft&i at ti C e)/2mR.b- After this step the 
tunneling along the perpendicular direction h/J\D ~ 0.5 s" 1 is negligible compared to 
the time needed for the experiment r exp = 0.17 s. The condition for being in the so called 
quasi two-dimensional regime which requires k&T, \x <C hum = Ti-2it -21 kHz — k-Q-1 fiK 
is also well fulfilled. 

The increased effective coupling g rs 2. 5 go leads to a slight expansion of the condensate 
in all directions, whereas the enhanced harmonic confinement u = l.2u) tends to com- 
press the sample. After all we fill about 40 discs, where the occupation of the central 
disc amounts to N 2 d ~ 4000. The chemical potential is fi = h ■ 2.7kHz = 133 nK. 
Subsequently the triangular 2D lattice is ramped up to the desired final lattice depth 
using a smooth sigmoidal ramp within 150 ms [36]. After a short hold time of 5 ms all 
optical potentials are switched off within less than 10 fis and the atomic cloud is imaged 
on a CCD camera after 21ms TOF. 

The lattice laser beams are derived from a frequency and intensity stabilized Ti:Sa laser, 
which delivers up to 1.3 W of power at a wavelength of 830 nm. Beam shaping telescopes 
focus the laser beams to Gaussian waists of 125 /im. Special care has to be devoted to 
the relative phases of the three laser beams generating the triangular lattice. Since 
optical fibers are used to deliver the optical power to the experiment fast servo loops 
eliminating any phase noise imprinted by acoustic Brillouin scattering inside the fibers 
are employed. This avoids heating caused by parametric excitation [71 137]. 
We calibrate the lattice depth by parametric excitation [37J from the first to the third 
band at sufficiently large laser power, where the bandwidth of the bands becomes negli- 
gibly small, and comparison of the observed frequencies with band structure calculations 
(Fig. [2b shows a full 2D band structure calculation). For this purpose we create ID 
lattices using two of the three laser beams and permute through all pairs of beams to 
create a potential as symmetric as possible. 
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The SF-MI transition of ultracold atoms confined in an optical lattice has been under 
vast investigation since first pioneering theoretical [381 EH] an d experimental [3DJ studies 
over the last decade. Many aspects and characteristic features of this quantum phase 
transition have been addressed such as the visibility of interference peaks [3 H l4"2l l4"3l 135] . 
the possibility of observing fine structure in the central peak [4"4"l |4"5] , probing of the 
energy gap in the MI phase [3D] and as a striking demonstration of strong correlations 
in the MI phase noise correlations have been observed in various experiments [46J. 
Moreover continuative experiments could directly map the wedding cake like shell 
structure of different single-site occupation numbers in the MI by high resolution 
tomography [37], high accuracy rf-spectroscopy [31] and direct in-situ imaging [49|. 
The effects of finite temperature [5DJ EH E2 E3] and dimensionality [53] have also been 
in the focus of interest. A detailed understanding of the MI state is an indispensable 
requirement for many fascinating experiments that have been proposed in the field of 
quantum information processing [55] and simulation of solid state systems [5]. 

Following the experimental sequence outlined above we have investigated the SF- 
MI transition in a triangular optical lattice for the first time. Fig [3^ shows a series of 
time-of-flight absorption images obtained by varying the lattice depth across the SF-MI 
transition. We have checked the reversibility of the ramp process by first bringing the 
system deep in the MI regime and subsequently ramping the lattice down to a depth 
corresponding to maximal visibility and holding the atoms for a short time comparable 
to the tunneling time J/h, after which we observe a very high recurrence of the visibility. 

4-1. Visibility 

The visibility, determined as described in [31] and averaged over 15 — 20 experimental 
realizations is shown in a (semi-logarithmic) plot versus the Hubbard parameter r\ in 
Fig. [3b. The visibility of the diffraction peaks corresponding to reciprocal lattice vec- 
tors start to smear out at a lattice depth of about 3.5 Er corresponding to a Hubbard 
parameter rj = U/zJ close to the expected value of rj = 5.8 [56] and drops to very 
small values for lattice depth Vq > 8E T . Note that the number of nearest neighbours z 
enters the calculation and is not consistent with any other number than z — 6, which is 
expected for the triangular lattice. In an inhomogeneous system a typical feature of a 
perfect atomic Mott-insulator is a wedding cake like structure of the occupation number 
per lattice site [39J. The relatively strong inhomogeneity in our experiment resulting 
from the large trapping frequencies of the dipole trap, suggests a rather pronounced 
structure of alternating MI and superfluid shells [4Tl |4~2| Wf\ l4~8l [57] with a central site 
occupation n of 3 — 4 atoms. Consequently reproducible kinks in the visibility are ob- 
served as also reported in [31] . The locations of the kinks fit well to the critical values of 
the Hubbard parameter 77™ for the transition to a MI with occupation number n derived 
from a mean- field model [56]. 
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Figure 3: (a) Absorption images of atoms released from a triangular optical lattice 
for various lattice depths given in units of E r . (b) Visibility of the interference peaks 
(blue circles) across the phase transition against U/zJ. Note that the visibility has been 
normalized to its maximum value of V max = 0.48. The visibility starts to drop around a 
lattice depth of 3.5 — 4 E r . Reproducible kinks marked by blue arrows may indicate the 
formation of Mott-shells. The condensate fraction (red triangles) decreases as the lattice 
is ramped to deeper values. Pronounced kinks coinciding with those observed for the 
visibility are observed. The condensate fraction vanishes more quickly than the visibility 
(Nq/N = indicated by red arrow). The FWHM of the central peak of the momentum 
distribution (yellow stars) is constant for low lattice depth. Around the phase transition 
it quickly starts to grow - also exhibiting small but reproducible kinks. The width of the 
incoherent background is plotted for completeness (green squares). Constant below the 
phase transition where it is a measure for the temperature (indicated by a green arrow), 
it starts to increase linearly with respect to U/zJ in the MI regime as expected from 
perturbation theory. Dashed vertical lines indicate the critical values for the transition 
to a Mott insulating state with n = 1, 2 and 3 respectively 



4-2. Condensate fraction 

As the visibility starts to drop the central as well as the diffraction peaks begin to 
broaden significantly as can also be deduced from FigSb. We have fitted a bimodal 
distribution with a fixed condensate width to the peripheral peaks of the interference 
pattern. In order to account for interaction induced broadening, data from the first 



| The condensate width has been determined using fits to data corresponding to low lattice depths 
well below the phase transition. 
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Brillouin zone has been excluded [531 EB]- After sufficiently long time-of-flight the 
density distribution consists of sharp interference peaks [35] caused by the condensed 
atoms which are assumed to exhibit the same parabolic shape as the external harmonic 
confinement. The other part of the bimodal density distribution is produced by thermal 
atoms whose distribution has been modeled employing a Gaussian. For small lattice 
depths this seems reasonable for a bimodal BEC in a harmonic trap perturbed by a 
periodic potential. For deep lattices the envelope of the density distribution is given by 
the Fourier transform of the Wannier function |u>(q)| 2 which is again well approximated 
by a Gaussian. The condensate fraction determined in this way is also plotted and 
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Figure 4: Bimodal fits to the peak structure of the absorption images. Data from 
the first Brillouin zone has been excluded to minimize interaction effects during the 
expansion. The data points consequently represent an average of all 6 satellite peaks 
(blue shaded area in inset). Note that every absorption picture has been averaged over 
20-30 experimental runs. From left to right the lattice depth increases, reaching the 
critical point for singly occupied lattice sites in the last picture. 



accurately follows the decrease of the visibility in shape for shallow to intermediate 
lattices. This is expected since the contrast of the interference pattern is intrinsically 
connected to the superfluid density [59] and therefore the condensate fraction. 

4-3. Residual visibility - Particle-hole excitations 

For lattice depths close to the expected SF-MI phase transition the superfluid fraction 
disappears while a small residual visibility persists even far in the MI regime. This resid- 
ual visibility in the MI regime has been attributed to excitations in terms of particle-hole 
pairs [HI H2| which should in principle be present at any finite hopping J. We have 
determined the coefficients a by fitting the model density derived from a first order 
perturbed MI state [HI [60] according to (n k ) = iV|w;(k)| 2 (l + £i a cos^k/k,)). The 
expected value a = 4J/Un(n + 1) fits surprisingly well the observed data for an average 
occupation number of n = 2 in spite of the rich shell structure that is expected in our 
experiment with occupation numbers ranging from 1 to 4 (see Figj5]). The error of the 
applied fit starts to increase significantly (indicated by the intersection of the two red 
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lines in Fig. Ob) at a lattice depth corresponding to the appearance of a finite superfluid 
fraction supporting the idea of particle-hole pairs as the origin of the residual visibility. 
Although recent simulations [52] show that ultracold thermal atoms also exhibit a 
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Figure 5: (a) Absorption images of ultracold atoms released from a triangular optical 
lattice for selected lattice depth above the MI transition point are shown together with 
fits of first order perturbation theory with only one free parameter a. (b) Experimentally 
obtained values for a together with the theoretical expected curve for an average 
occupation of (h) = 2. Note that a is proportional to 8J/U. The fit error Aa exhibits 
a distinct increase at a Hubbard parameter corresponding to the appearance of a finite 
superfluid fraction as indicated by the red arrow in Fig.0 The red solid lines are guides 
to the eye. 



significant interference pattern, our experimental parameters suggest that loading the 
atoms into intermediate deep lattices should even result in adiabatic cooling [61] . How- 
ever, the vanishing of the interference pattern has to be evaluated with great care for 
atomic samples with a non- negligible thermal fraction (see [6lJ[62j[63],[35]). By ramping 
back from a deep MI state to zero lattice depth and analysis of the density distribution 
after TOF we have checked that only negligible irreversible heating occurred for short 
to intermediate hold times. 



Note that the broadening of the diffraction peaks as well as the condensate fraction 
also show a kink like behavior supporting the observation derived from the visibility 
measurement [4"5| . 
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5. Spin dynamics in a triangular lattice 

We have performed measurements of spin dynamics in lattices of different depth and at 
different magnetic fields in order to gain insight into the physics of magnetic systems in 
periodic potentials. In principle one would expect a crossover from mean-field regime 
spin dynamics as described in e.g. |25j to few particle spin oscillations isolated at 
individual lattice sites [6] . It has been found for the latter case that coherent oscillations 
between two two-particle states ipi = \rriF — 0, 0) ipf — | — 1, 1) occur that can be 
readily explained within the framework of a Rabi-like model [6] . For atoms with F = 1 
in an optical lattice the effective Rabi frequency is composed out of a coupling strength 
Qif oc gi J \ w(r) | 4 dr given by the spin-dependent interaction between the two states and 
a detuning 5 = 5 + S(B 2 ). 5 is given by a part proportional to the quadratic Zeeman 
shift S(B 2 ) and an offset term determined by the difference in interaction energy of the 
initial and final state 5q cx gi/2. 

Here we focus on the investigation of spin dynamics at intermediate lattice depths which 
is associated with an overall spin independent trapping potential from the beginning. 
We present measurements of spin dynamics of 87 Rb BEC in F — 1 starting in an almost 
isotropic crossed optical dipole trap and ranging across the superfluid regime and up to 
the SF-MI phase transition inside a triangular optical lattice. In addition the magnetic 
field has been varied in order to change the relative influence of quadratic Zeeman effect 
and spin dependent interaction. We show that this ratio determines the behavior of 
oscillation amplitude and period. After a short reminder on the physics and notations 
of mean-field spinor physics within the single mode approximation and a description of 
the experimental sequence, we will present results suggesting a new and elegant method 
to tune the spin dynamics resonance known from experiments in harmonic traps (see 



5.1. Mean- field theory of spin dynamics 

The description of spinor condensates presented here is based on the evolution due to 
the corresponding mean-field Hamiltonian using the single mode approximation (SMA) 
as outlined in more detail in e.g. [251 126] . Spin mixing occurs as a consequence of 
two-particle collisions and can be classified by the total spin / of the colliding pair. For 
F = 1 spinor condensates the allowed values are / = 0, 2. The spin mixing interaction 
can be parameterized by the coupling coefficient gi = ^- a2 ~ a ° , each total spin channel 
being associated with an s-wave scattering length aj. The dynamics of the 2F + 1 
components of the vectorial spinor order parameter directly follows from the resulting 
mean-field energy functional for the spin part of the wave function (: 



e.g. m 128]). 



<#spin> = |iV(n)( 5l (F) 2 ) + hN(-p(F z ) + q(FD) 



(7) 



with 




(8) 
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Eff(W)«0 (9) 

ijfc 

|p| = |^-i?o and \q\ = parametrize the linear and quadratic Zeeman effect 
respectively. The dynamics governed by the above Hamiltonian can be solved exactly 
for F — 1 and a particular initial state C-/2 which can be easily prepared experimentally 
by the application of an initial n/2 rf-pulse j25j EE]. The solution of the resulting 
equation of motion for the populations of the individual m^-states can be expressed in 
terms of Jacobi elliptical functions (JEF) 

|Co(*)| 2 =(l-ksnl(qt))/2, (10) 

|C ± i(t)| 2 = (l + fcsn!(gt))/4. (11) 

The occurrence of the parameter k = gi(n)/q in the analytic solution indicates that 
amplitude and oscillation period of spin dynamics are determined by the ratio of 
interaction to quadratic Zeeman energy. The JEF in Equ. [TT] can further be simplified 
by approximations in two limiting cases: 

Zeeman regime (k <C 1) In the Zeeman regime, where the quadratic Zeeman effect 
dominates over spin-dependent interaction, the JEF can be approximated by 
ordinary trigonometric functions snfc(x) ~ sin(x), cnfc(x) ~ cos(x), dn^(x) ~ 1, 
yielding 

|Co(t)| 2 Z = R (l-A;sin 2 (gt))/2, (12) 

The amplitude of the oscillation is proportional to k, while the period is simply 
given by ir/q. 

Interaction regime (k ^> 1) In this regime the spin dependent interaction dominates 
any dynamics. The appropriate approximations for the JEF read 
snfc(x) = r sin 1 / fc (fcx) resulting in population oscillations 

|Co(*)| 2 '= (l-lAsin 2 Mn>t))/2. (13) 

Now the amplitude scales as 1/k and the spin-dependent interaction determines the 
period ir/gi(n). 

In the intermediate regime where both energies are of the same order of magnitude 
a spin dynamics resonance will occur which has been observed for the first time in a 
harmonically trapped 87 Rb F = 2 spinor condensate in [26J. 

[jj Fig.[6]shows measurements for different values of k for such a spin dynamics resonance. 
Note that the tuning of the resonance has been achieved by changing the magnetic field 
value while keeping the atomic density as constant as possible. 

§ F = 2 spinor condensates have an additional spin-dependent interaction parametrized by gi ■ |«Sb| 2 . 
However, due to the smallness of this contribution the resulting dynamics can well be approximated by 
the solutions found for F = 1. 
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Figure 6: (a) Spin dynamics resonance phenomenon in F = 2 taken from [26]. Plotted 
are coherent oscillations for different values of q/(g 1 (n)), where q ~ B 2 is the quadratic 
Zeeman energy and g\ (n) the first spin-dependent mean-field energy (processes involving 
Amp = 1). In this figure, gi(n) = 47s -1 (corresponding to a density of 1 x 10 14 cm~ 3 ) 
is used as a reference value, since only the magnetic field has been varied among the 
different data sets. The mean density has been obtained from SMA theory fits. 
Amplitude (b) and period (c) have been extracted for qualitative comparison with 
theoretical predictions for F = 1. 

5.2. Spin dynamics of S7 Rb F = 1 spinor condensates in a triangular lattice 

In this paper we discuss for the first time, that the mean-field spin dynamics resonance 
persists even in the presence of a periodic potential as long as a significant superfluid 
fraction is still present in the system. Moreover the renormalization of the spin- 
dependent interaction by the lattice potential allows for an effective tuning of the 
resonance at constant magnetic field. 

The experimental procedure for the results presented here is as follows. After 
preparation of BEC in |1,— 1) with particle numbers of N m 10 5 and no discernible 
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thermal fraction, we load the atoms in the three-dimensional optical lattice consisting 
of a 2D triangular lattice and a perpendicular ID standing wave as discussed above. For 
these measurements the individual lattice depths have been adapted to yield equal U/zJ 
for the different spatial directions. After the final lattice depth has been reached, we 
allow the system to relax to equilibrium for 5 ms before we apply a 7r/2-pulse to prepare 
the initial state Cn/2- The magnetic field is then abruptly switched to the desired final 
value. We then hold the atoms for a variable evolution time to allow spin dynamics to 
take place. In order to have atomic samples which are not optically thick after time-of- 
flight, we quickly ramp down the lattice within 5 ms following the end of the evolution 
time to obtain atomic clouds that allow for a reliable and unambiguous atom number 
determination. Finally all trapping potentials are switched off and after 21ms of time- 
of-flight including Stern- Gerlach separation the atoms are imaged with resonant light. 
In the superfluid regime the addition of a periodic potential should mainly show up 
in terms of an increased effective spin-dependent coupling owing to an increased mean 
density of the sample. We expect that the resulting effective k will determine frequency 
and amplitude of spin dynamics within this model. However different additional effects 
e.g. caused by reduced mobility or an increased thermalization rate might lead to 
deviations from this behavior and the data has been carefully analyzed if any indications 
for those deviations can be found. 

To estimate atomic densities in the different lattice regimes we have employed a Thomas- 
Fermi-like approximation of the Bose-Hubbard model for the superfluid regime and 
the regime close to the SF-MI phase transition. The corresponding densities are 
calculated by taking into account the particular overall harmonic trapping frequencies 
uj, onsite interaction matrix elements U and total particle numbers N according to 
n = 4/7 n = \ij g with 



Having obtained the mean density in this way a specific effective k value can be assigned 
to each set of measurements. It should be mentioned, that our experimental parameters 
allow for a tuning of the effective density by a factor of approximately 5 by changing 
the lattice depth. 

Fig. shows spin dynamics measurements in a crossed dipole trap (Hubbard parameter 
r\ — 0), in the superfluid regime (rj = 1.45) and in the vicinity of the SF-MI phase 
transition (rj = 8.9) for exemplary magnetic fields. It can be inferred from the graphs 
that the population in the m F = component p 00 undergoes periodic oscillations. The 
dynamics of the other components follows according to /0_i_i = p\\ = 1/2(1 — poo)- 
Oscillation amplitude and period clearly change depending on the lattice depth and the 
magnetic field. The time evolution of the mp = population has been fitted with a 
function p 00 = (1 — Asin 2 (B ■ (t + C)))/2 to extract the oscillation amplitude A and 
period 2tt/B. 

A direct comparison of amplitudes and periods obtained in this way to the analytical 
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Figure 7: (a) Exemplary measurements of spin mixing dynamics for different strength of 
the periodic potential and different magnetic fields. From top to bottom measurements 
have been performed in the dipole trap at \B\ = 270 mG, in the superfluid regime at 
\B\ = 30mG and close to the phase transition from superfluid to Mott-insulator at 
\B\ = 800 mG. The individual graphs show the time evolution of the population in 
the mp = state poo- Data is presented along with fits (solid lines) to determine 
oscillation amplitude and period, (b) Period of spin mixing dynamics as a function of 
the dimensionless parameter k c g = gi(n e s)/q for experiments in a crossed dipole trap 
(blue circles), in the superfluid regime in the presence of a triangular optical lattice (red 
triangles) and close to the quantum phase transition from superfluid to Mott-insulator 
(yellow stars). Solid lines are guide to the eye. (c) Theoretically expected period and 
amplitude according to Equ.[TTl 



mean-field solution Equ{TT] shows indeed good agreement if the spin-dependent 
interaction energy entering k = g± (n) /q is renormalized by the effective density in the 
optical lattice as explained above (Fig. [7b). A major consequence of this observation is 
that by raising an optical lattice potential, the atomic density and thereby the effective 
k can be varied which can readily be used to precisely tune the spin dynamics resonance 
at constant magnetic field and constant particle number. Beyond that we want to 
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emphasize that the physics of spinor condensates in optical lattices demands for further 
investigation. For example first indication for spin demixing at very low magnetic 
fields as well as its suppression with increasing lattice depth has been observed in our 
experiments. This effect leads to a slow-down and potentially stop of spin mixing 
dynamics. The accurate and detailed analysis of this and similar effects is however 
beyond the scope of this paper and wil be presented elsewhere. 

6. Conclusion and Outlook 

In conclusion we have presented first experiments with ultracold atoms loaded in a 
triangular optical lattice. In particular the MI transition has been investigated in a two 
dimensional system. In direct comparison to a cubic lattice the transition occurs at much 
lower lattice depth due to the much stronger localization of the atoms. Moreover we have 
observed spin dynamics of 87 Rb F — 1 spinor condensates in an optical lattice and found 
qualitative agreement with the mean-field solution raised in [25] as long as the lattice 
is shallow enough to maintain a significant superfluid fraction. The differences in the 
results at different lattice depths can mainly be reduced to an effective spin-dependent 
interaction gi(n aoTm ) which allows for a tuning of the spin dynamics resonance. 
Future perspectives include loading and investigation of spinor-BEC in the hexagonal 
anti-ferromagnetic lattice and analysis of magnetic ground states. Fascinatingly the sign 
of the tunneling matrix element J can be reversed by slightly modulating the frequencies 
of the three laser beams creating the triangular lattice as described in [20]. Thereby 
access to investigate fascinating model systems known from solid state physics that 
give rise to effects like Neel ordering and spin liquid phases is possible. Experiments 
concerning transport and spin-dependent tunneling for the creation of highly entangled 
multi-particle states are highly interesting research directions. 

7. Acknowledgements 

We thank the DFG for funding within the Forschergruppe FOR801 and Graduiertenkol- 
leg GrK1355. 

[1] I. Bloch, J. Dalibard, and W. Zwerger. Many-Body Physics with Ultracold Gases. Rev. Mod. 
Phys., 80:885, 2008. 

[2] M. Albicz, R. Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M. K. Oberthaler. Direct 

Observation of Tunneling and Nonlinear Self-Trapping in a Single Bosonic Josephson Junction. 

Phys. Rev. Lett, 95:010402, 2005. 
[3] A. Sorensen, L.-M. Duan, J. I. Cirac, and P. Zoller. Many-Particle Entanglement with Bose- 

Einstein Condensates. Nature, 409:63, Jan 2001. 
[4] C. Ospelkaus, S. Ospelkaus, L. Humbert, P. Ernst, K. Sengstock, and K. Bongs. Ultracold 

Heteronuclear Molecules in a 3D Optical Lattice. Phys. Rev. Lett, 97(12):120402, 2006. 
[5] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen, and U. Sen. Ultracold atomic gases 

in optical lattices: mimicking condensed matter physics and beyond. Advances in Physics, 

56:243-379, 2007. 

[6] A. Widera, F. Gerbier, S. Foiling, T. Gericke, O. Mandel, and I. Bloch. Coherent collisional spin 
dynamics in optical lattices. Phys. Rev. Lett., 95:190405, 2005. 



Ultracold quantum gases in triangular optical lattices 



17 



K. I. Petsas, A. B. Coates, and G. Grynberg. Crystallography of optical lattices. Phys. Rev. A, 
50, 1994. 

S. Trotzky, P. Cheinct, S. Foiling. M. Fcld, U. Schnorrbcrgcr, A. M. Rey, A. Polkonikov, E. A. 

Dcmlcr, M. D. Lukin, and I. Bloch. Timc-Rcsolvcd Observation and Control of Superexchange 

Interaction with Ultracold Atoms in Optical Lattices. Science, 319:295, 2008. 
P. Cheinet, S. Trotzky, M. Fcld, U. Schnorrbcrger, M. Moreno-Cardoner, S. Foiling, and I. Bloch. 

Counting Atoms Using Interaction Blockade in an Optical Lattice. Phys. Rev. Lett., 101:090404, 

2008. 

L. K. Lim, C. M. Smith, and A. Hemmerich. Staggered- Vortex Supcrffuid of Ultracold Bosons in 

a 3d Optical Lattice. Phys. Rev. Lett., 100:130402, 2008. 
Y.-J. Lin, R. L. Compton, A. R. Perry, W. D. Phillips, J. V. Porto, and I. B. Spielman. Bose- 

Einstein Condensate in a Uniform Light-Induced Vector Potential. Phys. Rev. Lett., 102:130401, 

2009. 

R. G. Melko, A. Paramckanti, A. A. Burkov, A. Vishwanath, D. N. Sheng, and L. Barents. 
Supersolid Order from Disorder: Hard-Core Bosons on the Triangular Lattices. Phys. Rev. 
Lett, 95:127207, 2005. 

L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, H. Fehrmann, and M. Lewenstein. Atomic 
Quantum Gases in Kagome Lattices. Phys. Rev. Lett., 93:030601, 2004. 

S. Wessel and M. Troyer. Supersolid Hard-Core Bosons on the Triangular Lattice. Phys. Rev. 
Lett, 95:127205, 2005. 

C. Wu, W. V. Liu, J. Moore, and S. Das Sarma. Quantum Stripe Ordering in Optical Lattices. 
Phys. Rev. Lett., 97:190406, 2006. 

L. Mathey, S.-W. Tsai, and A. H. Neto. Exotic superconducting phases of ulracold atom mixtures 
on triangular lattices. Phys. Rev. B, 75:174516, 2007. 

G. Juzcliunas, J. Ruseckas, M. Lindberg, L. Santos, and P. Ohbcrg. Quasirelativistic behaviour 
of cold atoms in light fields. Phys. Rev. A, 77:011802(R), 2008. 

L. H. Haddad and L. D. Carr. The nonlinear Dirac equation in Bose-Einstein condensates: 

Foundation and symmetries. Physica D, 238:1413-1421, 2009. 
S. Nakatsuji, Y. Nambu, H. Tonomura, O. Sakai, S. Jonas, C. Broholm, H. Tsunetsugu, Y. Qiu, 

and Y. Maeno. Spin Disorder on a Triangular Lattice. Science, 309:1697-1700, 2005. 
A. Eckard, P. Hauke, P. Soltan-Panahi, C. Becker, K. Sengstock, and M. Lewenstein. Frustrated 

quantum antiferromagnetism with ultracold bosons in a triangular lattice. arXiv, page 

0907.0423, 2009. 

H. -J. Miesner, D. M. Stamper-Kurn, J. Stenger, S. Inouye, A. P. Chikkatur, and W. Ketterle. 
Observation of metastable states in spinor Bose-Einstein condensates. Phys. Rev. Lett, 82:2228, 
1999. 

D. M. Stamper-Kurn, H.-J. Miesner, A. P. Chikkatur, S. Inouye, J. Stenger, and W. Ketterle. 
Quantum Tunneling across Spin Domains in a Bose-Einstein Condensate. Phys. Rev. Lett., 
83(4):661, 1999. 

M. Erhard, H. Schmaljohann, J. Kronjager, K. Bongs, and K. Sengstock. Bose-Einstein 
condensation at constant temperature. Phys. Rev. A, 70:031602(R), 2004. 

H. Schmaljohann, M. Erhard, J. Kronjager, M. Kottke, S. van Staa, L. Cacciapuoti, J. J. Arlt, 
K. Bongs, and K. Sengstock. Dynamics of F = 2 Spinor Bose-Einstein Condensates. Phys. Rev. 
Lett, 92:040402, 2004. 

J. Kronjager, C. Becker, M. Brinkmann, R. Walser, P. Navez, K. Bongs, and K. Sengstock. 
Evolution of a spinor condensate: coherent dynamics, dephasing and revivals. Phys. Rev. A, 
72:063619, 2005. 

J. Kronjager, C. Becker, P. Navez, K. Bongs, and K. Sengstock. Magnetically Tuned Spin 

Dynamics Resonance. Phys. Rev. Lett, 97:110404, 2006. 
M.-S. Chang, C. D. Hamleyand, M. D. Barrett, J. A. Sauer, K. M. Fortier, W. Zhang, L. You, and 

M. S. Chapman. Observation of Spinor Dynamics in Optically Trapped 87Rb Bose-Einstein 



Ultracold quantum gases in triangular optical lattices 18 

Condensates. Phys. Rev. Lett, 92:140403, 2004. 
[28] M.-S. Chang, Q. Qin, W. Zhang, L. You, and M. S. Chapman. Coherent spinor dynamics in a 

spin-1 Bose condensate. Nature Physics, 1:111, 2005. 
[29] W. Zhang, D. L. Zhou, M.-S. Chang, M.-S. Chapman, and L.You. Coherent spin mixing dynamics 

in a spin-1 atomic condensate. Phys. Rev. A, 72:013602, 2005. 
[30] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, and D. M. Stampcr-Kurn. Spontaneous 

symmetry breaking in a quenched ferromagnetic spinor Bosc-Einstcin condensate. Nature, 

443:312, 2006. 

[31] J. Kronjager, C. Becker, P. Soltan-Panahi, K. Bongs, and K. Sengstock. Spontaneous pattern 

formation in an anti-ferromagnetic quantum gas. arXiv, page 0904:2339, 2009. 
[32] W. Zhang, D. L. Zhou, M.-S. Chang, M.-S. Chapman, and L.You. Dynamical instability and 

domain formation in a spin-1 Bose-Einstein condensate. Phys. Rev. Lett., 95:180403, 2005. 
[33] M. Vengalattore, S. R. Laslie, J. Guzman, and D. M. Stamper-Kurn. Spontaneously Modulated 

Spin Textures in a Dipolar Spinor Bose-Einstein Condensate. Phys. Rev. Lett., 100, 2008. 
[34] F. Gerbier, A. Widera, S. Foiling, O. Mandel, and I. Bloch. Resonant control of spin dynamics in 

ultracold quantum gases by microwave dressing. Phys. Rev. A, 73:041602, 2006. 
[35] F. Gerbier, S. Trotzky, S. Foiling, U. Schnorrberger, J. D. Thompson, A. Widera, I. Bloch, L. Pollet, 

M. Troycr, B. Capogrosso-Sansone, N. V. Prokof'ev, and B. V. Svistunov. Expansion of a 

Quantum Gas Released from an Optical Lattice. Phys. Rev. Lett., 101(15):155303, 2008. 
[36] D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams, and P. Zoller. Creation of a Molecular Condensate 

by Dynamically Melting a Mott Insulator. Phys. Rev. Lett, 89:040202, 2002. 
[37] T. A. Savard, K. M. O'Hara, and J. E. Thomas. Laser-noise-induced heating in far-off resonant 

optical traps. Phys. Rev. A, 56:R1095-1098, 1997. 
[38] M. P. A. Fisher, P. B. Weichman, G. Grinstcin, and D. S. Fisher. Boson localization and the 

superfluid-insulator transition. Phys. Rev. B, 40:546-570, 1989. 
[39] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller. Cold Bosonic Atoms in Optical 

Lattices. Phys. Rev. Lett., 81:3108-3111, 1998. 
[40] M. Greiner, O. Mandel, T. Esslingcr, T. W. Hansen, and I. Bloch. Quantum phase transition from 

a superfluid to a Mott insulator in a gas of ultracold atoms. Nature, 415:39-44, 2002. 
[41] F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, and I. Bloch. Phase Coherence of an 

Atomic Mott Insulator. Phys. Rev. Lett, 95:050404, 2005. 
[42] F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, and I. Bloch. Interfernce pattern and 

visibility of a Mott insulator. Phys. Rev. A, 72:053606, 2005. 
[43] W. Zwerger. Mott-Hubbard transition of cold atoms in optical lattices. J. Opt. B, 5:S9-S16, 2003. 
[44] V. A. Kashurnikov, N. V. Prokof'ev, and B. V. Svistunov. Revealing the superfhiid-Mott-insulator 

transition in an optical lattice. Phys. Rev. A, 66:031601, 2002. 
[45] S. Wessels, F. Alet, M. Troyer, and G. G. Batrouni. Quantum monte carlo simulations of confined 

bosonic atoms in optical lattices. Phys. Rev. A, 70:053615, 2004. 
[46] S. Foiling, F. Gerbier, A. Widera, O. Mandel, T. Gericke, and I. Bloch. Spatial quantum noise 

interferometry in expanding ultracold atom clouds. Nature, 434:481-484, 2005. 
[47] S. Foiling, A. Widera, T. Miiller, F. Gerbier, and I. Bloch. Formation of Spatial Shell Structure 

in the Superfluid to Mott Insulator Transition. Phys. Rev. Lett, 97:060403, 2006. 
[48] G. K. Campbell, J. Mun, M. Boyd, P. Medley, A. E. Lconhardt, L. G. Marcassa, D. E. Pritchard, 

and W. Ketterle. Imaging the Mott Insulator Shells by Using Atomic Clock Shifts. Science, 

313:649-652, 2006. 

[49] N. Gcmelke, X. Zhang, C. L. Hung, and C. Chin. In situ observation of incompressible Mott- 

insulating domains in ultracold atomic gases. Nature, 460:995, 2009. 
[50] F. Gerbier. Boson Mott Insulators at Finite Temperatures. Phys. Rev. Lett, 99:120405, 2007. 
[51] G. Pupillo, C. J. Williams, and N. V. Prokof'ev. Effects of finite temperature on the Mott-insulator 

state. Phys. Rev. A, 73:013408, 2006. 
[52] Y. Kato, Q. Zhou, N. Kawashima, and N. Trivedi. Sharp peaks in the momentum distribution of 



Ultracold quantum gases in triangular optical lattices 



19 



bosons in optical lattices in the normal state. Nature Physics, 4:617-621, 2008. 
W. Yi, G.-D. Lin, and L.-M. Duan. Signal of Bose-Einstein condensation in an optical lattice at 

finite temperature. Phys. Rev. A, 76:031602, 2007. 
T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslingcr. Transition from a Strongly Interacting 

ID Superfluid to a Mott Insulator. Phys. Rev. Lett, 92:130403, 2004. 
T. Calarco, U. Dorncr, P. S. Julienne, C. J. Williams, and P. Zollcr. Quantum computations with 

atoms in optical lattices: Marker qubits and molecular interactions. Phys. Rev. A, 70:012306, 

2004. 

D. van Oostcn, P. van der Straten, and H. T. C. Stoof. Quantum phases in an optical lattice. 

Phys. Rev. A, 63:053601, 2001. 
K. Mitra, C. J. Williams, and C. A. R. Sa de Mclo. Superfluid and mott-insulating shells of bosons 

in harmonically confined optical lattices. Phys. Rev. A, 77:033607, 2008. 
I. B. Spiclman, W. D. Phillips, and J. V. Porto. Condenstate Fraction in a 2D Bose Gas Measured 

across the Mott-Insulator Transition. Phys. Rev. Lett, 100:120402, 2008. 
R. Roth and K. Burnett. Superfluidity and interference pattern of ultracold atoms in optical 

lattices. Phys. Rev. A, 67:031602, 2003. 
I. B. Spiclman, W. D. Phillips, and J. V. Porto. Mott-insulator Transition in a Two-Dimensional 

Atomic Bose gas. Phys. Rev. Lett, 98:080404, 2007. 
P. B. Blakie and J. V. Porto. Adiabatic loading of bosons into optical lattices. Phys. Rev. A, 

69:013603, 2004. 

B. DcMarco, C. Lannert, S. Vishveshwara, and T.-C Wei. Structure and stability of mott-insulator 
shells of bosons trapped in an optical lattice. Phys. Rev. A, 71:063601, 2005. 

T.-L. Ho and Q. Zhou. Intrinsic Heating and Cooling in Adiabatic Processes for Bosons in Optical 
Lattices. Phys. Rev. Lett, 99:120404, 2007. 



